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We argue that for generic systems close to a critical point, an extended Fluctuation-Dissipation 
relation connects the low frequency non-linear (cubic) susceptibility to the four-point correlation 
OO ' function. In glassy systems, the latter contains interesting information on the heterogeneity and co- 

operativity of the dynamics. Our result suggests that if the abrupt slowing down of glassy materials 
is indeed accompanied by the growth of a cooperative length I, then the non-linear, 2>uj response to 
an oscillating field (at frequency uj) should substantially increase and give direct information on the 
temperature (or density) dependence of I. The analysis of the non-linear compressibility or the di- 
^ . electric susceptibility in supercooled liquids, or the non-linear magnetic susceptibility in spin-glasses, 

should give access to a cooperative length scale, that grows as the temperature is decreased or as 
the age of the system increases. Our theoretical analysis holds exactly within the Mode-Coupling 
. Theory of glasses. 
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I. INTRODUCTION 



A yet unexplained property of fragile glasses is the extremely fast rise of their relaxation time (or viscosity) as 
the temperature is lowered, much faster than predicted by a simple thermal activation formula 0. If interpreted 
^ in terms of an effective activation energy, the latter increases by a factor five to ten between 1.5 Tg and the glass 
transition temperature Tg. The basic mechanism for this increase is not well understood, but it is reasonable to think 
^ ' that it is intimately related to cooperative effects 0, and possibly to the presence of an underlying critical point 
^ ! [slI^M-IMIgl lTollTlllT^ I. The dynamics becomes sluggish and the activation energy increases because larger and larger 
QQ I regions of the material have to move in a correlated way to allow for a substantial motion of individual particles. Long 
\^ ' time scales must be somehow associated with large length scales. Although the idea of a cooperative length has been 
\C) . discussed in the context of glasses for many years 0, 13 ' is only recently that proper measures of cooperativity 



(and of the size of the rearranging regions) were proposed theoretically (see fo'^ earlier insights) and measured 
in numerical simulations [T^ lisl ll^ ll^ (see also |3 'l^; 'is] for related experimental work) . The idea is to measure 
how the dynamics is correlated in space; technically, this involves a four-point correlation function which measures 
the spatial correlations of the temporal correlation (see Eq. ^ below for a more precise definition) . Recent extensive 
"j^ numerical evaluations of this four-point correlation function in Lennard- Jones systems have confirmed the existence 
of a growing length scale as temperature is decreased 12, _1£, JJj, and have shown that different observables, such 
I ■ as the relaxation time or the diffusion constant, scale as powers of this length, emphasizing its crucial importance 
'■^ ■ as far as the physics is concerned. In the framework of granular systems, diverging length scales near the jamming 
^ , transition have also been reported in numerical studies of model systems |20l. 

O • Although many different theoretical approaches to the glass transition IS Hi S ^i ITTl IT^ can potentially 
explain the existence of such a growing dynamical correlation length, these theories lead to rather different quanti- 
tative predictions for the behaviour of the four-point correlation function (see UJ). Thus, experiments measuring 
directly this four-point function would be extremely valuable to refine our understanding of the glass phenomenon 
^ and prune down the number of candidate models. Up to now, unfortunately, only indirect exp erimental indications 
' of a cooperative length scale associated to heterogeneous dynamics have been reported 0, Il8l |. 

On a different front, that of spin-glasses, length scale ideas have also been expressed in the recent years to account 
for non equilibrium phenomena such as aging, memory and rejuvenation effects |22, 23, 24, 25, 26, 27]. Although 
spin-glass order is not easy to define nor to detect, the idea is that some kind of domain growth occurs, whereby spin- 
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glass correlations establish on larger and larger length scales as the age of the system increases. The growth of this 
"coherence length" has been established numerically by comparing two replicas of the same system I29ll30ll3lll3^ . 
This trick is obviously inaccessible to experimentalists, who have nevertheless provided indirect evidence of a growing 
length scale, and some indications on its rate of growth with time and temperature psl Is^ 133. 135^. Again, a direct 
measure of this length scale is lacking - finding a clear-cut experimental signal of a cooperative length in disordered, 
amorphous systems would certainly be a major breakthrough [T3 |. 

The aim of this paper is to point out that in slow glassy systems at equilibrium, the non-linear (cubic) response 
to an external field (electric, magnetic, pressure, etc.) in fact probes directly the four-point correlation function 
mentioned above, and therefore the cooperative length it may contain. Our main prediction, detailed below, is that 
the 3lj harmonic response to an a.c. field of frequency uj and amplitude h is given by Xsi^j T) h^, where the non-linear 
susceptibility xs behaves at low frequency as: 

^,{^u,T)^^e-m{uT). (1) 

In the above relation, Xs is the static linear susceptibility, TL a certain complex function that depends weakly on 
temperature, and r is the temperature dependent relaxation time of the system, which can be directly measured 
using the linear susceptibility. The cooperative length £ (measured in units of the microscopic length ^ obtained 
from the point correlation function) is expected to grow as the temperature is reduced, and rj an exponent related to 
the spatial structure of the four-point correlation function |67j| . The above prediction holds for equilibrium systems; 
we will furthermore see below that in the case of glasses and spin-glasses in a field, Ti{{)) = 0. Below the glass 
transition temperature, on the other hand, the system by definition falls out of equilibrium. Its dynamics becomes 
non stationary and exhibits aging, which means that the effective relaxation time of the system increases with the 
age tw of the system [36, 37]. This increase of the relaxation time is again most probably related to the growth of a 
coherence length in the system, = £{tw)- Assuming simple aging behaviour, the generalization of the equilibrium 
result CQ) then reads: 

X,iu;,t^) = ^el-^Tiojt^), (2) 

which should allow one to extract from non-linear aging susceptibilities a non equilibrium coherence length, in a 
much more direct way than previous attempts. (In the above equation, J- is another scaling function, which also 
contains possible violations of the standard Fluctuation Dissipation Theorem and the appearance of a non trivial, 
dependent, effective temperature ^3)- central results, Eqs. (|1I2(I . that we will motivate below, states that (a) 
the non-linear susceptibility has the same frequency scaling as the linear susceptibility, which is not surprising and (b) 
it grows as the cooperative length increases, which should allow a direct experimental test of the relationship between 
length and time scales in glassy systems. 

As for comparison with previous works, the divergence of the static non-linear susceptibility at the spin-glass (at 
zero field) or dipolar-glass transition, displayed by Eq. at w = 0, is of course well documented, both theoretically 
|42| and experimen tally .42. .43,. .44] . The generalization to the dynamical non- linear susceptibility in the critical region 
was also discussed |23, |43, Eg but not, to the best of our knowledge, its generalization to the non-equilibrium, aging 
regime, Eq. ((SJ. The situation for glass- forming liquids is quite different, since no static phase transition with a 
diverging static susceptibility has ever been identified, neither in experiments nor in simulations. Purely based on an 
analogy with spin-glasses, it was suggested in j47| that the non-linear dielectric constant of molecular glasses rnight 
grow as the glass phase is approached (although this was not borne out by the experiments done at that time [43 )■ 
A similar suggestion was made in concerning the non-linear compressibility of soft sphere binary mixtures, with 
numerical results that are not incompatible with a substantial increase of X3i^ ~ 0) ^ the temperature is lowered. 
We will show below that such a growth is indeed expected, although the theoretical situation for glass formers is much 
less clear than for spin-glasses - in particular, xs although is growing may never diverge in glass formers. Different 
scenarii for the glass transition can be envisaged and lead to quite different predictions, for example on the value of 
ry and on the relationship between i and r (or tw ) ■ 

In the following section we will give some physical arguments that motivate our results, and muster the predictions 
of different theoretical models for glass-formers. A more detailed and technical discussion is then presented in Section 
III. Finally our conclusions are presented in Section IV. 
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II. PHYSICAL ARGUMENTS AND RESULTS 



A. Spin-glasses 



1. Order parameter and non-linear susceptibility 



Let us first focus on spin-glasses in zero external magnetic field, H = 0. These systems are known, both theoretically 
and experimentally, to have a non zero transition temperature below which the magnetization profile, (s^), freezes into 
one (or more) amorphous configurations. The ordered state is characterized by a non zero Edwards- Anderson (EA) 
parameter q = [(sx)'^], where (...) indicates thermal averaging and the brackets a spatial (or disorder) average. These 
systems display an unusual type of long-range order, which cannot be detected using either one body or two-body 
spin-spin correlations: because the ordering pattern is random, the average magnetization [{sx}] remains zero and the 
spin-spin correlation [{sxSy)] short-ranged, even in the spin-glass phase. Correspondingly, the linear susceptibility, 
related by a fluctuation dissipation theorem (fdt) to the integral of the spin-spin correlation function, does not diverge 
as Tg is approached, even if some long-ranged correlations appear in the system. The way to get rid of the spurious 
cancellation between strongly correlated and strongly anti-correlated spins is well known: exactly as one should square 
(sx) to obtain a non zero Edwards- Anderson parameter, one should also square (sxSy) before averaging over disorder. 
The integral over space of that quantity now diverges as Tg is approached, and in fact has two interesting physical 
interpretations. The first one is the susceptibility of the spin-glass order parameter to small random ordering fields. 
Imagine one adds small random magnetic fields on every site. Using linear response, one can write, for a given 
sample at T > Tg and H = 0: 



where the subscript means that the correlation functions are evaluated at zero external field. Squaring this relation, 
summing over x and averaging over the random fields gives the sensitivity of the EA order parameter to a random 
pinning field: 



Clearly, the divergence of XSG signals an incipient instability towards an ordering pattern favoured by the small 
pinning fields, exactly as the divergence of the usual two-body susceptibility signals an instability to ferromagnetic 
order, triggered by a small (uniform) magnetic field. 

As defined above, xsG has a clear theoretical interpretation but seems hard to access experimentally. Fortunately, 
there is a direct relation between xsG and the non-linear susceptibility, which can be directly measured. The intuitive 
idea is that the non-linear susceptibility is actually a measure of the (quadratic) dependence of the linear susceptibility 
on the external field. Using fdt the change of the connected correlation function between two spins (and hence of 
the linear susceptibility) induced by the field contains the term: 



Averaging over space (or over disorder), only the terms {z = y, x, z' = a;, y) survive, the first one giving \(sxSy)^ as in 
XSG- A more precise treatment for Ising spins at zero field leads to the exact relation X'i{'^ = 0) = ^(3xsG~2)/(A;Br)'^. 
Therefore, the static non-linear susceptibility of spin-glasses diverges at the spin-glass transition temperature, a well- 
known effect that allows one to measure some of the critical exponents experimentally [42. i43j . The physics behind 
the correlation induced amplification of X3 is clear: the influence of the polarization of spin Sx on Sy may be either 
positive or negative, but it has the same sign as the reverse influence of Sy on Sx- Therefore, the quadratic effect of 
an external field h on the dynamical correlation between any pair of spins has a well defined sign, in turn leading to a 
diverging non-linear susceptibility as the size of correlated regions increases, even if the linear susceptibility remains 
small. 




(3) 



y 




(4) 




(5) 



2. Non zero external field 



The case where a non zero external field H is present is more subtle. In mean-field, the spin-glass phase survives 
in a whole region of the T, H plane, below the de Almeida-Thouless (at) line. The spin-spin correlation function 



4 



[{{sxSy) — {sx){sy))'^] is long-ranged in the whole spin-glass phase but is no longer directly related to the static non- 
linear susceptibility. Some exact compensation mechanism |49l l52 | actually cancels the divergence in the combination 
of four-spin correlations appearing in Xsi'^ — 0)- Therefore, the non- linear susceptibility is finite in the whole spin- 
glass phase. There is in particular no divergence of Xsi^ = 0) on the AT line, except at H — 0; rather, the non-linear 
susceptibility is discontinuous across the AT phase transition j5Qj. Within the droplet theory, on the other hand, the 
spin-glass is detroyed by any non zero field; both the spin-glass and non-linear susceptibilities are finite when H =/= 
|22l l5lj. For H = 0, a, compensation mechanism similar to that of mean-field glasses is also at play, but does not 
prevent the non- linear susceptibility to diverge for all T <Tg |2^ . 



3. Dynamical non-linear susceptibility 



The above qualitative arguments for the static non-linear susceptibility can be extended to the dynamical case as 
well. As will be recalled below, the dynamical fdt gives: 

1 f d 

^j^Y.j dh — {s.Ati)sy{t:i))ohy[h). (6) 

Therefore, the change in the connected dynamical correlations between Sx{ti) and Sy{t2), induced by a uniform, but 
time dependent external field, will contain a term like: 



/d d 
dt3dt4 — {Sx{tl)s^{t3))o — {Sy{t2)s^>{t4))oh{t3)h{t4 
dt-i dU 



(7) 



Repeating the same argument developed in the static case, i.e. averaging over space (or disorder) and using fdt to 
relate the connected correlation function to the dynamical linear susceptibility leads to a non-linear response function 
that reads: 

^3 

X3{h;t2,h,U) - ^ rft dt dt [^^^(^i)^i/fe))o(s:rft4)syft2))o]- (8) 
y ^ ^ 

Taking ^1,^2,^3,^4 all within an interval of the order of the relaxation time t of the system, we see that the correlation 
function entering ^3 above defines a cooperative length scale £, such that the dynamics of Sx and Sy within this time 
interval is dominated by common events. This in turn leads to our scaling prediction, Eq. (Q, near the transition 
temperature. The exact result for the dynamical X3(ti; t2, ^3, ^4) needs to be worked out carefully (see Section 111), 
since fdt for higher order correlations is more involved than for two point functions |5,l| . Although different from 
the above naive expression, it indeed contains four-spin correlation functions that capture the cooperativity of the 
dynamics. Intuitively, again, the non-linear response is strong if on the scale of the relaxation time, two spins feedback 
on each other's dynamics - this cross correlation is squared and survives averaging, even if the correlation itself is of 
random sign. 



4- The aging regime 



In the low temperature, spin-glass phase, the relaxation time r is infinite and the age of the system tyj plays a 
crucial role - all time dependent correlation functions depend explicitly on tw |37l |. However, exciting the system 
with a field of frequency 1 /tw will give the non- linear response of a spin-glass equilibrated only up to a certain length 
scale = £{t^) |63|- Interestingly, contrarily to standard ferromagnets, spin-glasses are thought to be critical in 
their whole low temperature phase, in the sense that the space integral of the connected correlation fimction [{sxSy)'^] 
diverges for all T < Tg. Within the mean- field replica theory, the static non- linear susceptibility within one phase is, 
as mentioned above, finite (except for T — Tg and H = 0) We however think that the mechanism cancelling 

the divergence does not operate at finite frequencies, and that an equation similar Eq. |(2J) will hold in the aging 
phase, but with an infinite number of time domains rather than the simple scaling variable Lut^ 36, 372- The explicit 
calculation of in the context of a spherical p-spin model would be extremely interesting; in particular one may ask 
whether the effective temperature appearing in the non-linear response is the same as that appearing in the linear 
response 55] . 

In the droplet picture [S^, |^, on the other hand, the static non- linear susceptibility diverges for all T < Tg, 
provided H = 0, and one should certainly observe a non-linear susceptibility increasing as in Eq. ||2Jl, although the 
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original droplet model with activated dynamics would rather predict a function T of Inw/lnfm and a logarithmic 
growth of l(t^~). The peak value X'i^w-,^ = l/^iu) should grow as The numerical value of 2 — 77 is yet 

unknown, but following Fisher and Huse [23 1 one may expect d — 36 < 2 — fJ < d — 9 with 9 « 0.2 in three dimensions. 

B. Structural glasses 

1. Four-point density functions 

Let us now discuss structural glasses. The important lesson we learn from spin-glasses is that a non-trivial amor- 
phous type of long-range order can set in. In the case of glasses, the subtlety comes from the absence of quenched 
disorder; however, there has now been many papers exploring the idea of self-induced disorder which could drive a 
similar transition in homogeneous, frustrated systems (see, e.g. [37| and refs. therein). This has led, in particular, 
to the "Random First Order Transition" scenario (5|, where a glass transition of the same nature as the spin-glass 
transition in mean-field p-spin models takes place (see for recent quantitative results) . Whether a true transition 
of this type can exist in real systems with finite range interactions is still an actively debated issue; it is nevertheless 
extremely fruitful to explore the properties of systems for which this transition is, in some sense, nearly present. The 
order parameter in the would-be glass phase is the amplitude of the frozen in (random) density fluctuations Sp,j: . As 
for spin-glasses, the average of this quantity is zero, but {Spx)"^ is not, and plays the role of the Edwards- Anderson 
parameter. Similarly, one expects (SpxSpy) not to show any interesting features (beyond that typical of a liquid 
structure factor), whereas its square may reveal long-range cooperative dynamics. The analog of the spin-glass and 
non-linear susceptibilities discussed previously can be easily found in the case of glasses: the former can be seen as 
the susceptibility to a random external pinning field js^ that triggers the appearance of one particular type of frozen 
density fluctuation, whereas the latter is directly related to the non-linear compressibility, i.e. the response to a 
pressure field that couples to the density. Other non-linear responses to a field that couples to the degrees of freedom 
undergoing a glass transition are also relevant (for example, the dielectric response when the dipoles are strongly 
coupled to the translational degrees of freedom, such as in glycerol, OTP, etc.). 

Let us directly focus on the dynamical susceptibility and postpone the discussion of its static limit to section 
II. B. 3. Indeed, the analogy with spin-glasses has to be taken with a grain of salt (see section II. B. 3). The dynamical 
four-point density function defined as: 

G4(r, t) = {5px{t = 0)5px{t)5px+r{t = 0)(5p,+,(i)) - C^{t) C{t) = (<5p,(t = 0)(5p,(t)), (9) 

is related to the dynamical non-linear response of the system to an external excitation that couples to the density 
|69|| . Once again, the idea is that the change of the two point correlation between x and y induced by the external 
field of frequency uo ^ t^^ will be large if on that time scale, the dynamics at these two points is strongly correlated, 
which is true precisely if Gi{x — y, t) is large. The extended, non- linear fdt discussed in the next section makes this 
statement more precise and finally leads to our central results, Eqs. H1I26|I ,7Qj . Now, recent numerical |,1S..15...16, 1'^ 
and theoretical work j^, flOl l2ll l58| have focused on the above choice of four point density function. Its integral 
over space X4(t) = / d'^rGi{r,t) (divided by V) gives the variance of the correlation function C{t) for a system of 
finite volume V [56l |. and is therefore a good quantitative measure for dynamical heterogeneities. This quantity was 
unambiguously shown to display a peak at t = r, of increasing amplitude as the temperature is decreased and the 
glass temperature is approached [ij, [la, Ha. Il7| , signaling increased cooperativity in the dynamics and the growth of 
a length scale t, which should in turn show up in the non-linear response of the system. 

2. Different scenarii for the glass transition: qualitative predictions 

We therefore expect, on very general grounds, the non-linear response to a field that couples to degrees of freedom 
undergoing a collective freezing phenomenon, to increase substantially (as ^^~'') as the glass phase is approached. 
However, as we discuss now, the details of this increase do depend on the specific scenario at play. Most important 
in that respect is the quantitative relation between the cooperative length scale i and the relaxation time r, which is 
often a power-law t ^ where z is the dynamical exponent. 

One scenario for the glass state is based on the idea that some mobility defects are needed to trigger the dynamics, 
which slows down at low temperatures because these defects become rare ^ 58, 59] . Kinetically constrained models 
provide an interesting framework to quantify this idea. The class of so-called "East" models seems able to capture 
some of the phenomenology of fragile glasses and predict a temperature dependent exponent z = zq/T^ which may 
become large at low temperature [59l |60| . This is very important since the relaxation time of the system is known to 
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increase by 15 orders of magnitude as the temperature is reduced from 1.5Tg to Tg. But if, say, z = 15 the cooperative 
length £ would only increase by a very moderate factor 10. 

Another well known scenario is based on the Mode-Coupling theory (mct) of supercooled liquids, which predicts a 
dynamical singularity at a finite temperature Tmct, where the relaxation time should diverge as r ~ (T — Tmct)~^ 
|61i |. with a non-universal exponent 7. It was recently shown that this singularity is actually accompanied by the 
divergence of a cooperative length precisely defined in terms of the four-density correlation function above, Eq. 
The nature of the transition is actually equivalent to that of the mean- field p— spin glass, where both the 

spin-glass and the non-linear susceptibility diverge at T^j^j, [Tlj . The exponent z is found to be equal to 27; 
reasonable values of z being in the range 4 — 6. However, the ideal MCT phase transition is 'avoided' in real systems. 
Only the first 2 — 4 decades of increase of r can be satisfactorily accounted by mct, before some new physics come into 
play, that smear out the mct transition. In the temperature region where T < Tmct^ the system should, according 
to mean-field, be completely frozen. In finite dimensions, however, barriers to motion are finite and the dynamics is 
instead strongly activated. More precisely, the liquid is a 'mosaic' of local metastable glass states, that must rearrange 
collectively 00, [nl. The size of these frozen clusters is the cooperative length ^, which increases as T decreases, but 
now only logarithmically with r. Therefore, within the Random First Order Theory (rfot) of |3| which unifies mct 
and the mosaic scenario one expects a cross-over between a power-law increase of £ for T ~ Tmct and a much more 
modest increase of £ as the temperature is reduced from Tmct to Tg. [Below Tg, aging effects come into play and we 
expect that an equation similar to Eq. Q will hold in that regime.] Finally, the 'avoided critical point' scenario of 
Kivelson and Tarjus also predicts a cooperative length that grows weakly (logarithmically) with the relaxation time 

. . . 

The value of 77 in Eq. above is not known either, and presumably depends both on the scenario and on the 

temperature regime. Negative values of 77 « —1.58 have been reported for East models 59] , whereas rj is probably small 
in the Mode-Coupling region |lOj] . In the simplest mosaic state scenario where clusters are compact, the exponent 2 — rJ 
is equal to the dimension of space d. Rather independently of the scenario, we therefore expect a noticeable rise of the 
non-linear susceptibility in supercooled liquids as the temperature is reduced: taking ^ = and £(Tg)/£(1.5Tg) — 5 
leads to an increase of the peak value of Xd. by a factor « 25. (Note however that depending on the particular physical 
observable there might be other contributions coming from the temperature dependence of Xs or, for example for the 
non-linear dielectric susceptibility, from the Lorentz field effects that may play an important role in strongly dielectric 
systems [6^'). 

3. Glasses vs. spin- glasses: some caveats 

The tricky aspect of the analogy between glasses and spin-glasses is that the static non-linear susceptibility of glasses 
is in fact not expected to display any divergence. As a matter of fact, no growing correlation length has ever yet been 
found in any static correlation function close to the glass transition, neither in experiments nor in simulations. This 
is not only because the glass transition is never sharp in real glasses, but also because from a purely theoretical point 
of view the situation is more subtle. Within rfot, for example, the static non-linear susceptibility does not diverge. 
The reason is that rfot predicts an exponentially large number of possible amorphous states and thus a non zero 
configurational entropy for Tk <T < Tmct |3,f3^, "4^ (where Tr- is the 'entropy crisis' temperature). Now, since 
equilibrium thermal averages are sums over all states, one has {SpxSpy) = J2a'^a{Spx^Py)a, where the subscript a 
indicates that the average is restricted to the metastable state a, and Wa is the weight of state a. Therefore: 

^{SpxSpy)^ = ^ ^ WaWi3{SpxSpy)a{SpxSpy)i3. (10) 
y a, 13 y 

The divergence of the static spin-glass susceptibility is due to the diagonal terms a = f3. In the case of spin-glasses, 
the number of relevant states is effectively finite, and the above sum diverges in the spin-glass phase. Within rfot, 
on the other hand, the number of relevant metastable states is so huge for Tr- < T < Tmct that the diagonal 
contribution tends to zero, as does the off-diagonal contribution since the different frozen patterns are uncorrelated 
with each other. The only way to unveil any growing correlation is to focus on the terms a = /3, a calculation that 
is possible in mean- field (see previous footnote). In non mean- field situations, this can be done in two ways: one can 
study (a) the static correlations but in a self-induced static pinning field (freezing all the particles outside a cavity 
and studying the thermodynamics inside |llj]'). which select a particular state a or (b) the dynamical correlations 
on a time scale short enough for the system to remain in a single state 0. For example, the four-point correlation 
function defined in Eq. ^ for t ~ r, effectively reproduces, for r large enough, the static sum restricted to a = /3 |^, 
since r is the time needed to evolve from state a to any another state. More precisely, one sees that in this regime. 
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6*4 can be written in a form closer to the corresponding expression for spin-glasses: 

where the disorder average [...] in the case of spin-glasses is replaced by an average over time, [...]t', corresponding 
to different configurations of the self-induced disorder in glasses. Therefore, we expect that the analogy with static 
spin-glasses indeed makes sense for t ^ t. A practical consequence of this observation is that the scaling function 
TC{x) in Eq. should tend to zero at zero frequency for glasses, at variance with the spin-glass case where it remains 
finite at T = r+ and for zero field. 

Even the spin-glass case turns out to be tricky since, as mentioned above, the scaling function Ti.{x — > 0) in fact 
also goes to zero for H ^ in the context of the full replica symmetry breaking solution [i^ Is^ . Nevertheless, our 
dynamical result Eq. (|25|l suggests this sum rule will generically not hold at finite frequency whenever the four-point 
correlation function has a non trivial time dependence. This statement should of course be checked explicitely for 
mean-field models with continuous replica symmetry breaking [H^. If indeed Ti{x ~ 1) is found to be non zero close to 
the AT line, the experimental study of the dynamical non linear susceptibility, predicted to diverge for lot ~ 1, would 
offer a direct way to prove or disprove the existence of an AT line in real systems (see for a recent discussion) . 



III. THEORETICAL ANALYSIS 



In the rest of this paper, we give some theoretical justifications of our central result, Eqs. H1I2I) . We will use 
the Langevin equation formalism for continuous spins j^, but our result are expected to hold more generally (for 
example, if the continuous spins are replaced by interacting particles with Newtonian dynamics). 



A. Linear response 



We assume that the equation of motion of spin Si is given by: 

dtS,^-d,,H + r],{t), (12) 

where H is the Hamiltonian of the system, which we do not specify explicitly. In the case of spin glasses it contains 
quenched disorder and possibly one body terms ensuring an Ising like character to the spins s^. The coupling to an 
external, site dependent field hi{t) amounts to add to H the sum over spins of hi{t)si. The Gaussian noise rji is as 
usual of zero mean, white in time and decorrelated from spin to spin: 

{miti)rijit2)} = 2kBTSih - t2)6,,j. (13) 

Since the noise is Gaussian, one can establish the following identity: 

,/ dsdh) 



^dhHm-ksT^^^). (14) 



Let us first quickly re-establish the standard linear fdt. From the above equation and the equation of motion, the 
response of a spin to an earlier field is: 

X^jituh) - (1^) = j^{^dhm.,s,+dsMih)]) . (15) 

The averaging above assumes the system to be in equilibrium: we average over all histories with initial conditions 
appearing with the equilibrium Boltzmann weight. The second term in the right hand side is zero since, for an 
arbitrary observable 0{{ta}) that depends on times ta, all posterior to t2, one has: 

{0{{U})ds^H{t2)) = -keT fYlds{ta)P[{sita)}\s2]Oi{sita)})ds2ds,e,^p[-H{s2)/kBT] = 0, (16) 

a 

where S2 = 5(^2) and the last equality holds because the last term is a total derivative. Therefore, one finds the usual 
FDT relation: 

x»iti,t2) = ;fc^^(«*(^i)^^(^2)). (17) 
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Integrating this quantity over t2 with a constant field hi(t2) = h gives the static susceptibihty Xs, which, as is well 
known, is found to be the integral over space of the two-body correlation function. In the case of a static critical 
point where the correlation length ^ diverges, one would have Xs ~ ^^^^"''/^b^; where 77 is the standard critical 
exponent of the static transition and fi the elementary magnetic moment. However, in the case of glassy systems, the 
two point function is not critical and one rather expects Xs ^ IJ-^'^/ksT where ^ remains microscopic and does not 
grow appreciably lowering the temperature (or increasing the density). As emphasized in Section II, one should in 
the case of amorphous systems rather focus on non linear effects to observe some non trivial behaviour. 



B. Non-linear response: the static limit 



As a consequence we want to extend the above calculation to the response at time ti > t2 to three field 'kicks' at 
times ^2 > ^3 > ^4- This is given by: 

X3,..(^l,t„t3,t,) ^ ( ,,^.(,^g:(tUft4) ) = (^-^)-^(^^(*^)^.(*^)^^(^3)..(i4)). (18) 

Using three times the Langevin equation of motion, and once the above trick to get rid of the final dsiH{t4), we find 
the following general relation, involving four terms: 

{kBTfx3.ijkl{tl,t2,t3,t4) = {Si{ti)Sj{t2)sk{t3)si{t4)) + 

CLt2Cit5Ctl4 

S^{tl)^s^Hit2)sk{t3)slit4)) + —— {s^itl)sJ {t2)ds, H {h)si{t4)) + 



dt'idt4 ^ dt2dt4 

±{s,{t4)ds^H{t2)ds,H{h)si{t4)). (19) 

Let us first analyze the static limit of this expression. From the above result, one can show in full generality that the 
static non- linear susceptibility xss = Xii^ — 0: T), obtained by integrating over all ^2 > > ^4 with a constant field 
hiit) = ft, on all sites, is given by: 

{kBTfx^s - ^ ^(s.(ti)s,(ti)sfc(ti)si(ii))c, (20) 

ijkl 

where the subscript c means that one takes the connected part of the correlation and N the total number of sites. 
This result is exact and can be obtained directly using equilibrium statistical mechanics. In the present context, only 
the first term in expression Eq. H19|l for X3,ijfc; (^i: ^2, ^3, ^4) contributes for a; = 0. As discussed in Section II, the 
long range order setting in spin-glasses is unveiled not by the two-body correlation that oscillates in sign and averages 
to zero, but by the square of this two-body correlation. Therefore, the leading dominant term in the above sum 
corresponds to the square of the two-body correlation obtained pairing k, I - say - i with j and k with I within 
the two-body correlation length ^ (which typically remains small at all temperatures): 

(fcsr)3x3s«-^e''E^'fc G^k^{s^{tl)sk{t,))l. (21) 

ik 

If one now assumes that Gik scales as in usual critical phenomena [2^ l49j | : 

' gf^V (22) 



then the sum over i, k behaves as N£'^^'^, finally leading to a static non-linear susceptibility given by 



X3s 



^^"^ - 7^^^"" (23) 



(fcBT)3 ksT 



where C is a constant, and i is counted in units of the static correlation length ^. This is the zero frequency result 
given in Eq. Note that for spin-glasses in a non zero field, a full replica symmetry breaking calculation reveals 
that Gik = {si{ti)sk(ti))^ is short-ranged which means that Xs is in fact non divergent. 
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C. Non-linear response: dynamics 

The extension to non zero frequency of the above result can proceed in different ways. Our resuh Eq. Q can be 
sim ply seen as a standard dynamical scaling assumption close to a critical point, as is indeed correct for spin-glasses 
I45I l46l | . This result is expected to hold whenever a critical point is responsible for the simultaneous increase of the 
relaxation time and the cooperative length. This is true of the Mode-Coupling Theory of glasses 0,0], and also of 
other scenarii discussed in the introduction and in Section II, which rely on the existence of an underlying critical 
point [S IE Si E| ■ From a more technical point of view we want to justify that the behaviour of the non-linear 
cubic response is the same as of the first term on the right hand side of Eq. H19|l , whereas the three other terms are 
either negligible or of the same order of magnitude (on frequencies of the order of r^^), but not more divergent. 

A simple case that can be treated in some generality is when the fluctuation of the norm of the spins can be 
neglected, for example for Ising spins that can be recovered from the Langevin equation in the limit of infinitely 
sharp double well potential that is zero if = 1 and infinite otherwise. After several integration by parts and using 

= 1, one can show that the three last terms of Eq. (|19|l do not contribute to the non-linear a. c. susceptibility 
at low frequencies (much smaller than the microscopic, high frequency scale of the model). One is therefore left with 
the first term of Eq. (|19|1 . that contains three derivatives with respect to time. If one assumes that the four-body 
correlation {si{ti)si{t2)sk(t3)skit4)) is, for \i ~ k\ £, only a function of (ti — ^2)/^, (^2 — ^3)/''' a-nd (^3 — tA)/T, 
the integration over t2,t3,t4 with an oscillating field at frequency u) and over space directly leads to Eq. i.e. a 
non-linear susceptibility that scales as a certain function Ti of ujt. This result is only justified in the low frequency 
domain; for high frequency, contributions from the short-time /3-regime will obviously come into play. Note that very 
generally, we expect Ti to be non trivial, although it does vanish at zero frequency whenever the static susceptibility 
is finite, as is the case for glasses and spin-glasses in an external field (see the discussion in II. A. 2, II. B. 3). 

More generally, one can argue both physically and diagrammatically that the three last terms of Eq. H19|) give 
contributions which are at most of the same order of magnitude as the first one. From a physical point of view, these 
terms contain less time derivatives that the first, but also contain the local 'force' acting on the configuration, dsH{t). 
Since we are interested in the low frequency response of the system, we can decompose the dynamics of the spins 
into a fast part and a slow part s*, that corresponds to the dynamics on a time of order r. It is clear that the 
force acting on the slow modes can only lead to a slow dynamics of these modes, i.e. \ds-'H\ ^ t^^ . Therefore, for 
frequencies ~ r~^, one has, for example, 

{s^{.h)dsMt2W{h)si{U)) ^ T-'F (^^1^, ^i^) , (24) 

(where F is a certain function) , which after integration leads again to a result of the form . 

One can understand this result from a different point of view using diagrams for a general Langevin equation, which 
leads to a dynamical field theory with the spin field s and the response field s 0| . The non- linear cubic response xs 
at time ti to three instantaneous fields at times ^2,^3,^4 can be written in full generality as (see Fig. 

X3iti;t2,h,U) = J <dt^dtj,<r,,,,,,,(t;;i'2,t^,,t^)x(ii -t'i)x(i2 -i2)x(i3 -^3)x(i4 -^4) (25) 
where F^^s^s^s is the amputated vertex with legs s, s, s, s (for simplicity we skip here the space indices). Note that the 
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FIG. 2; Diagrammatic representation of G4 



vertex rs^s,s,s is zero by causality because it contains for sure a closed loop of response functions. The other vertices 
do not appear because the correlation function (ss) vanishes by causality. 

Now let us consider the diagrams contributing to the connected four body correlation function. There is a first 
contribution Gf obtained by plugging three two body correlation functions into Ts^s,s,s, see Fig. 12 It is straightforward 
to check, using fdt, that this first series of diagrams, G4, is directly related to the non-linear cubic response. If Gf 
was the only contribution, then one would find an extended fdt where only the first term on the right hand side of 
Eq. (fT^ contributes. 



There is another contribution, G\, that corresponds to constructing ladders with the irre- 
ducible vertices ^Yl-s rs^^Ts-s s^^fs-s s^^fs-s s^^fs-s s left of F^.s^^^s and the irreducible vertices 
r^r;«.,,r-r., ,,F^, ;,%^, :,r^^, :,r-^.,_:,' on the 'right of r,,,,,,,. See fig! El [We recaU that the irre- 



ducible vertex TY2-3 4 is the sum of all Feynman diagrams contributing to Fi,2;3,4 (the amputated vertex) that has 
the property that cutting two internal lines does not separate the diagram into two disconnected parts, such that one 
part contains the lines 1, 2 and the other one the lines 3, 4.] Finally, the last contribution, G4, is formed by plugging 
together the same irreducible diagrams used in G\ but without making use of Ts^s,s,s- 

In the case of the critical equilibrium dynamics of spin-glasses, Eq. Q is already known [ZM li^ and could have 
been guessed a priori from the general scaling properties of second order phase transitions with a single diverging 
length (and time) scale. As a consequence in this case G\, G% are of the same order or less divergent than G4. The case 
of structural glasses is a priori more tricky, since there is no consensus on the effective critical microscopic model that 
would describe them. However, if we take as an established fact (at least numerically) that the four body correlation 
is governed by a length scale that increases as the glass is approached, then this effect has to be contained in (at least) 
one of the three contributions G4,G4,G4. Now, the non-linear response certainly contains the contribution related 
to G4; therefore both xa and G4 grow (or even diverge) similarly unless another family of (more) diverging diagrams 
(the ones contributing to G4,G4) can be constructed. We believe that this is a rather unlikely scenario and instead 
we expect that in general Gf and G4 are of the same order of magnitude, and G4 is sub-dominant, in which case Eq. 
(fT|l holds. Strictly speaking, these arguments prove that if the non-linear cubic dynamical response diverges, a similar 
(or stronger) divergence is expected for the four body correlation function, but not vice-versa. Therefore it would be 
important to check our prediction for specific models of the glass transition, in the spirit of [23|. Here, we just want 
to emphasize that the techniques used in can be used to establish that our central result, Eq. holds within 
the Mode-Coupling Theory of the glass transition. 

Finally, let us remark that the extension to the non-equilibrium case can be tackled in a similar way. In partic- 
ular, since the four-body correlation function diverges with ^S^] in spin-glasses, and the classification in terms of 
G4,G4,G4 carries over to the non-equilibrium case, the above discussion can be generalized to the non-equilibrium 
case as well. 

In summary, we have shown in this section that for glassy systems close to a critical point, where the relaxation time 
and cooperative length diverge, an extended approximate fdt relates the non-linear susceptibility to the four-point 
correlation function in the low frequency domain: 



dt2dt3dti 



{si{ti)sj{t2)sk{h)si{t4)); 



(26) 



where ~ means that right and left hand side have the same critical behaviour. The additional terms missing in the 
above equation are either of the same order of magnitude, or negligible. 
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FIG. 3: Diagrammatic representation of G4 
IV. CONCLUSION 

In conclusion, we have shown in this paper that if the abrupt slowing down of glassy materials is accompanied by 
the growth of a cooperative length then the non-linear, Zlo response to an oscillating field (at frequency uj) should 
substantially increase and give precious information on the temperature (or density) dependence of The theoretical 
motivation is that the non-linear susceptibility is approximatively related, for glassy systems close to a critical point, 
to the four-point correlation function that captures dynamical cooperativity. This relation is certainly correct within 
the context of the Mode-Coupling Theory of glasses, but should hold in other cases as well. 

In supercooled liquids, the analysis of the non-linear compressibility (sound wave harmonics) should allow one to 
probe directly the existence of a growing cooperative length. This should also be true of the non-linear dielectric 
susceptibility, at least in systems where the dipoles are strongly coupled to the glassy degrees of freedom. Although 
early experiments seemed to show no interesting effects [43, we believe that more systematic studies should be 
performed |62l |. especially now that numerical simulations have unambiguously shown the growth of a cooperative 
length in the four-point function [T^ITsLITslll^lr^ . These experiments should also allow one to bridge the gap between 
the length-scales observed on simulation time scales and the length-scales observed experimentally on much larger 
time-scales close to the glass transition temperature The study of non- linear specific heat effects, although more 
complex, may be interesting too |65l l66j | . From a more general point of view any non- linear dynamical response (for 
example, non-linear rheology in soft glassy materials) should be worth studying if the corresponding linear response 
can be used as a probe of slow dynamics. 

In spin-glasses, non-linear a.c. magnetic susceptibility measurements in non-zero field could shed light on the 
existence of a de Almeida-Thouless line. In the aging phase, such measurements should allow one to test in more 
details the length scale ideas put forward in [2^ I23 . T24L \Wl l26l . Compared to the case of glasses, the experimental 
situation is particularly encouraging since the non-linear susceptibility is already known to diverge at the spin-glass 
transition. There should be a clear trace of this divergence in the aging phase, except if some subtle cancellation 
occurs even at non zero frequency (the mean-field replica theory indeed predicts such a cancellation in the static 
case). The effect of temperature cycling on the non-linear susceptibility should then give direct indications of the 
mechanisms of rejuvenation and memory [25l l26l | . We therefore hope that the ideas expressed in this paper will help 
shed light on the issue of dynamical heterogeneity and cooperativity in disordered, amorphous systems. 
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